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ANNOTATION 


Wo present a numerical analysis of 1 nonlinear axisymmebric our 1 ** 
rents in a viscous liquid in spherical annulil of various thickness 
when only the inner sphere rotates* 

We obtained the lines of flow and the lines of equal angular 
velocity s the spatial spectra of the kinetic, energy and of individual 
velocity components j we calculated the integral characteristics of 
the motions 5 the total kinetic energy and the torque for various 
values of the numbers Re* 

In a thin annulus » the nonuniquenes 3 of the steady-state solu- 
tions of the Navier-Stokes equations was established* and the regions 
(with respect to Re) in which various flow regimes exist were 
determined* 

Calculations of the currents in a thick annulus together with 
the results of previous experiments make it possible to draw conclu- 
sions concerning the stability of currents in the Re range considered. 

Comparison of the numerical results with experiments shows good 
quantitative and qualitative agreement , 


Nonlinear Axisymmetrio Liquid Currents in 
Spherical Annuli i 

H. M, Astaf'yeva, N, D* Vvedenskaya, and 1. M, Yavorskaya 


INTRODUCTION 


The study of nonlinear shearing currents in a liquid or a gas in /5 » 
rotating spherical annulus is the basis for an understanding of many 
dynamic processes of global scale in oceans, the atmospheres of the 
Earth, planets, and in the interior of stars of various spectral 
classes. 

In all these instances, the form of the currents and their 
stability depends essentially on two important factors: the spherical 
geometry of the volume in which they flow, and the rotation. Neglect 
of even one of these factors can lead to results quite remote from 
reality, as experiments show [1 - 6]. 

One of the simplest shearing currents in which such factors are 
considered is a current in a viscous incompressible liquid in a 
spherical annulus, the boundaries of which rotate around the same 
axis with different constant angular velocities. By analogy with the 
plane and cylindrical cases, we call it the Couette spherical current . 

The published results of experiments [1 - 6] have shown that, 
depending on the values of the parameters of similar problems: 
the Reynolds number Re-Cn^r^v, the relative thickness of the 

s 
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* and the Rossby number Q, 


(or the 


annulus v * w — _ 

z i £d 

Reynolds number with respect to the external sphere Rg^r; ) * /6 




the character and form of the currents arising in the annulus turn 


out to be extraordinarily diverse (Zf? inland Sk are the radii 
and the angular velocities of rotation of the internal and external 
spheres* respectively* and \> is the kinematic viscosity). 


For sufficiently small numbers Re [7 - 10]* the Couette spherical 
current is axisymmetric and consists of differential rotations around 
the axis and the meridional circulation whose intensity and form for 
fixed numbers Re are determined by the parameters 6 and e. Here it 
is already clear that the analogy between the plane and cylindrical 
currents* on the one hand* and the spherical* on the other, is purely 
external: the first two currents are one-dimensional and have one 
velocity component in the direction of the displacement of the velo- 
city. A current in a spherical annulus is much more complex: it has 
three velocity components and depends on the two spherical coordi- 
nates r and 6* and on the Reynolds number. 

This difference becomes still more significant when the basic 
current loses stability. In the plane case* if the very first in- 
stability leads to the spontaneous development of turbulence and* 
in the cylindrical case* to the appearance of stationary axisymmetric 
or (when the Rossby numbers are large and negative) periodic undulat- 
ing Taylor vortices* then in the spherical annulus [5* 11] when the 

i 

basic current loses stability, there arise diverse currents which 
differ widely among themselves* depending on the thickness of the 
annulus and the number e. Obviously this indicates that there are 
various mechanisms which cause the instability of the basic current 
for different 6 and e. 


Such a difference between plane, cylindrical and spherical 

currents arouses not only applied* but also theoretical interest in 
« 

the study of Couette spherical currents* since they obviously are 
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more general in character. This greater generality leads, as is to 
be expected, to significantly greater mathematical difficulties in- 
obtaining a solution: indeed, for arbitrary values of similar 

parameters, calculation of the spherical current reduces to the 
solution of a boundary value problem for a system of nonlinear 
partial differential equations. 

It is possible, therefore, that the most interesting results 
were those obtained experimentally in recent years. The analytical 
and numerical results are few in number and relate only to certain 
specific points in the space of similar parameters, A detailed 
analysis of all results in papers published up to 1974 is given in 
our survey article [12]. Here we shall touch on only those results 
in recent years, of which our paper is a direct development. 

First of all, let us note that spherical annulii can be sub- 
divided into thin annulii where 6 <6*, and thick annulii with 
5 > 6 s , depending on whether or not the principle of “change of 
stability" holds at the limit of stability when He = Re , The criti- 

w 

cal value of the Reynolds number is a function of 5 and e, and the 
critical value of the thickness of the annulii 5*, for which the 
character of the secondary current changes also, is obviously a 
function of the Rossby number e. Henceforth, we shall consider only 
rotations of the internal sphere for which 6* = 0.24. 

Up the the present, the question as to the stability of the 
Couette spherical current in thick annilii has not been resolved 
unambiguously. For a long time it was assumed that thick annulii 
were always stable [1, 2, 13]* The fact which the energy theory [103 
of the qualitative difference between thick and thin shells did not 
discover is that the lower boundary of the limit of stability in 
thick annulii was somewhat higher than in thin annulii. According 
to this theory, the most dangerous disturbances, for all 6 and e 
turned out to be the nonaxis ymmetr leal disturbances with azimuthal 
wave number m = 1. 


■ i / / 

r 

V 

The most detailed experiments have shown, however, that thick /S 
annulii are also unstable „ In the papers [3, 6], the authors demon- 
strated visually' the loss of stability by currents in annulii of 

thickness 0*24 £ I for Re c s , significantly exceeding the 

critical values of He E from the energy theory and exceeding the anti- 
cipated critical values obtained by extrapolation from the experi- 
mental stability curve in thin annulii. At the limit of st ability * 
a complex secondary nonaxisymmetrical current was observed, con- 
sisting of the same number of vortices in each hemisphere and 
propagated In an- azimuthal direction at a certain phase velocity 

Ql($< Q-i • 

i 

A completely different result was obtain in [4], where the 

stability of the currents in thick annulii with S’ »X f 27 and 2,25 
were studied by measuring the torque M transmitted to the external 
sphere, and by visual means. On the curve M -M(Re), four breaks 

were found in the current in the annulus with S' =1,37 , the first" 

three when the numbers R.0 ^ were smaller than 

R,S C (S" s=I , 33} (obtained in [6, 11]), and only the fourth break 

occurred when *1,33). Usually, breaks in the integral 

characteristics of the currents correspond to the onset of Insta- 
bility and the emergence of a secondary current visually distinct 
from the pre-critical current. Visually, no significant rearrange- 
ments of the current were detected for the first three breaks; the 

fourth break, corresponding to 5*1,33) of [11], was Identi- 

fied with a sudden complete turbulization of the flow. In [4] It 

was proposed that for Rg^ the first break corresponds to a drastic 

instability of the basic current of finite amplitude, which leads 
to secondary currents of the same form as the basic current. 

Additional experiments conducted by us [11] in a thick annulus 
with s*X,33 with the help of the film thermal flowmeter ”DISA U /£ 


4 


mounted in the equatorial plane- of the current, have shown that 
there are no changes in the characteristics of the current up to the 
values * When the recorded signal, its spectral 

density and the correlation function have shown the presence of only 
one frequency of the visually observed vortices* Moreover, investi- 
gation of the spectral characteristics of the current for the highly 

supercritical numbers have also shown that there are 

in the current no more than two different oscillation frequencies. 
Thus it was' established that the transition through Re » Re c , corre- 
sponding to the fourth break, is not accompanied by the sudden emer- 
gence of turbulence, but by the transition to a complex stratified 
periodic motion such as may be found in [ 3 ]- The question regarding 
the first three breaks remained open, since the change to a 
steady-state regime established with the aid of a thermal flow- 
meter is much more complicated. 

\ 

To clarify this question in the present paper, numerical 
investigations of liquid currents were undertaken in an annulus of 

thickness § = 1,33 in the interval 0 & ££50, the results of 

which are disucssed in sections 2 and 3 . YJe obtained curves of the 
torque and the kinetic energy as functions of the Reynolds number, 
and we studied the structure and scale of the prevailing motions. 

Previous numerical investigations of the currents in a thick 
annulus (£► ss I) were presented in [7 - 10]. However, the form 
in which these results were presented and the small number of them 
made them unsuitable for our purposes. 

In thin annulii with 4 0,24 in the supercritical region, it 
was discovered experimentally that there exist several steady-state 
axis ymme trie current regimes [2, 4 - 6, l*i, 15 ] such that some of 
them can exist for the same number Re, depending on the prehistory 
of the process — which agrees with the nonuniqueness of the solu- 
tion of the steady-state boundary problem for the Navier-Stokes 
equations. In particular, it was established experimentally that 


/IQ 
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in an annulus with thickness 33 Q,X3C,, there exist 4 steady-state 
currents with two* four, six, and eight annular vortices, parallel 
to the equator [6]. 

In the present paper, we study numerically all the experimen- 
tally obtained nonlinear steady-state axisymmetric currents in an 

annulus of S » O t Xl ; we establish the limits of their existence 
and the possible transitions from one regime to another* In con- 
trast to the experiment, numerical modeling made it possible to 
study the structure of the currents, their energy characteristics, 
and the linear scale of the prevalent motions, etc* 

We know of only one paper [14] in which the authors discuss the 
results of separate calculations for the flow functions of the cur- 
rents when <5 « ~X in thin annulii with § ss 0,178 for Re - 700, 900, 

and . 1500 and with b = 0,9375 ff « 8000 *» 8000 without any 
analysis of the sequential transition^ of the various regimes and 
the regions of existence and nonuniqueness of the currents. 

We present the results of calculations of the currents in a 
liquid in a thin and a thick annulus during rotation of the internal 
sphere only; the physical interpretation of the results and a com- 
parison with experiment are discussed later in sections 2 ’and 3- 

1, STATEMENT OF THE PROBLEM AND 

THE METHOD OF SOLUTION 
* 

The Couette spherical current is described mathematically by 
the solution of the boundary value problem for the Navier-Sfcokes 
system of equations. We shall find the steady-state solution by the 
method of adjustment, i.e., as the limit (if it exists) when t -*■ 
of the solution of the initial-boundary value problem' for the non- 
steady-state equations (t is time). In dimensionless form, the /ii 

problem has the form: 


6 




U VJ.U ss ~ 

V 7 - 


Caa) 




Ut'» u 


i wn 

I , i 


^ •5‘S^h^ I t when f 


U a ^Oj. If ^ ® & (**^'*& |%? 


(a.a) 


.0 *' • 1 

7' =L/ ix - 

/%0 V J ® fe*0 


V 

^ 


(a. 3 ) 


here* u is the velocity vector* p is the pressure! W* ^ 

are the components of u in spherical coordinates <p. For the 

scales of length, time, velocity and pressure, we assume, respect 
tively, the radius of the internal sphere the reciprocal of the 
angular rotational velocity of the internal sphere 

dj'ti and £p 0 tl *e density of the liquid), 

CO s* £Sq /^ 2-t *j£. 

For the numerical solution of the problem, we shall follow the 
procedure in [16]. 

We seen an ax i symmetrical solution of (2,1) - (2.3)* Let us 
represent the desired functions in the form of a series in terms of 

the associated functions °(&0 and ^0^) * where X - Ct96 Q , 

and the coefficients depend on r and t. We shall restrict ourselves 
to a finite number of terms in the series: * 

f A . . i - 7 11 . 3 \ A / i 


Ujt&t ) ~<L U&tfflefah %&)&(*) 

,u . 

. 4«. I d J PjfchS £ 


(a. 4) 


? 



hero in s the term ^ which is isolated in the 

steady^state solution of (2*1) ** (2*2) when Rg^O (the stokes 
current)* For this tern we have: 







<v ' 

j giii .T L i^ai|..i.i(i jpwniij'IW rM** rrt * 

y (H-i-^^i 


hot us represent the nonlinear terms in analogous form [for 
brevity all calculations are hereafter presented much as the first 
of the equations (2*1)3: 







Q 
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Let us especially stress the choice of the method of approximating t 

the terms ^ » 4 ** X* 8* S* We shall consider our expansions as 
interpolationsala l*e* s we shall require that the equations (2*p> 


be satisfied for 






Such a choice of Qj is explained by the fact that it is 
natural to treat our functions like the function QC ir CO, ^ Q , 

£ X 4?^, and the points 5C * are here the zeros of a 

d : ■ « 

Chebyshev polynomial. It is well-known that interpolation with 
respect to such a system of nodes is optimal in some sense.* 
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Substituting (2.4) and (2.5) into (2.1), we may now separate 

the variables, and obtain for Ua } V<i^ the following system 

of equations: * *’ 4 

VUa -f \ ■ (9 ,T>U* Z QUe l -£(?**) < 

■Orb Re Yi x ZZ Q>t t <dZ Z 3 - 


ftp 

n>? z , 

Z z 

= > 

| 

qvi . 
nrb 

~Te | 

> j_ 2L _ 

[ r*- 9z “ vz | 


z 

i J 5 i 

uuniL 1 \ 

‘i,fl ’ \ 

We 

n>^ 

i 1 

Re t 

"V I 

i 

^ n 

e !<? 
<?< 

Ole- 


£*• 


r **■ 


( 2 . 6 ) 
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l‘he boundary conditions assume the forms 

U&fat ) = V e h,i)^ Wg hi)=0 for z-i'^ '[+$. (2.7). 

It is clear that equations (2.5) are exactly satisfied on the rays; 


& - &. 


a «*a 
(/ ~ 
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Let us now carry out the discretisation of the problem with 
respect to t and r and pass from (2.6) to finite difference equa- 
tions. As a preliminary step, let us bake a change in the inde*» 

pendent variable r by means of the formulas J-60 , so that, 
after obtaining the equally spaced points of the difference scheme, 
the calculated points could he obtained which are not equally 
spaced relative to r. 


Thus, we shall consider all the functions only at the moment 
Vlt' p and at the points 6^- ft) fcK , 

k = 0* 1*0 yc 3 or at the points ij^y ~j (<) + k ? ic .-0 * 

The functions Uz ? Ve } We will be given at the points 

'ire,, and at the points • Let g (*tvi >\') ~ ^ Q) * The 

finite difference equations have the form: 



Here, e [&) is the coordinate in the expansion for ,8yt) 




uV*nufo) 

■ Um) QQ 



The* boundary conditions for (2.3): 



t v 



everywhere. 


It should be emphasised that in ('2,8) the linear terras are com- 
puted on the (n + l) st annulus, and the nonlinear right side per- 
tains to the h ' ^ annulus, so that on each new annulus a linear 
boundary value problem is, solved* 


We shall consider only the current which' is symmetric with 

respect to the equator r Q' ~ l therefore, in the expansion (2.*J) 

is is sufficient to restrict ourselves to the terms of corresponding 

parity. Namely, 

4 • ■ - 


Z U e 9= e ° i P = ’ 


u«t — 21 Oq 3% j 
€*0,9. ; / I 
l-i 


l-X 


^T\"C 


2l. Y<s./^^ j M-v ~~ 


1/^Sg V^3r5? 


In the calculations, L Is odd. The equations must be satisfied on 

X raysi 1 , $ • . _ . >l-f- 


(The choice of even L is Inconvenient, since in this case the system 
(2.6) has an additional degeneracy for &~L). The values of h for 
various versions are specified In section 2. 

The choice of any even step with respect to t requires a certain 
accuracy ; since our scheme is “explicit-implicit*** the inertial terms t 
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in it are calculated on the lower annulus* The estirnate from above 
for v which guarantees the stability of the scheme is obtained 
wlthoutv’dlffloulty in case the coefficients in (2.85 do not depend 
on k* It has the form: 

«tr < a m<^ju l^/Re * 

■ . i 

whore 0 *v 1* In the case of variable coefficients, the character of 
the estimate is preserved, but varies somewhat. For values of x in 
our calculations* of. the supplement. 


While establishing the solutions with respect to time, we varied 
the torque transmitted to the- external sphere; 

f ■ ■ 



‘and the kinetic energy of the current; 

* • Ti ,U$ . , \ „ . 

£ - 3? J o 1 (U^+Uf+Uj) % N **0 #© 


This was necessary so that these integral characteristics may be /16 
maintained with a precision of up to 0.1$. 


2'. RESULTS OP THE NUMERICAL CALCULATIONS 

Calculations of the boundary value problem for the Navier« 

stakes equation in a thick annulus with % » J,88 * when only the 

internal sphere was rotated, were carried out in the interval 

R© » § ISO in five steps relative to Re 5 close to the values 

and steps relative to Re were reduced to two. In 

addition, the oux^rent was calculated for 'Re * 1^0, 2QQ, and 22Q. 

For the initial data, we usually took |the values of the functions 
for the preceding step relative to Re . An example of the c&laula^ 
felon for Re » 200 is given in Figure 
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• In all the intervals of the Re values considered, a unique 

steady-state axisymmetric current was obtained, To verify the 
uniqueness of the solution, calculations with different initial 
data were performed. Thus, for example when p£^ ^ 55 

disturbances were introduced into the basic current computed for 
L - 31, disturbances which substantially change the scale of the 
prevailing motions: 1) the amplitudes of all harmonics, beginning 

with & > 8, were significantly increased, 4 and for & < 8, they re- 
mained unchanged; 2) the amplitudes ojf the harmonics with & < 4 
were equated to zero, and with & > 6, they Increased significantly. 

In all cases, the solution was qJickly established for the basic 
regime. 


In proportion to the growth in the Reynolds number, the number 
of associated functions L, used during the calculations, increased. 
Simultaneously with the characteristics of the current, the spectra 
of the kinetic energy and the individual components of the Velocity 
were calculated, depending on the meridional wave number l. The 
spectra made it possible to determine the order of the neglected 
terms, and thus to evaluate the accuracy with which the solution is 
approximated by the finite series (2.4). 


Since for Re = 60, it turns out to be sufficient to take L « 17, /l? 
and for Re = 150 to set L = 25, so that the kinetic energy of the 
neglected terms does not exceed 0.1# of the total kinetic energy of 
the current. Simultaneously with the calculation of the structure 
of the current, the integral characteristics were calcualted: the 
torque transmitted to the external sphere M (Figure 2), and the 
total kinetic energy E as a function :of Re (Figure 3). 

In a thin annulus, numerous calculations of the current were 
carried out from Re = 300 co Re < 2600. The concrete values of Re 

for which solutions were obtained are marked with dashes on the dia- 
gram (Figure 4). Usually when making the calculations, we took for 
the initial data the values of the current chax*aeterist±cs for the 
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preceding step relative to Re, .the step relative to Re being 
chosen positive as well as negative. 

One of the basic questions in the case of numerical calcula- 
tions by the direct as well as the semi-direct methods is the ques- 
tion concerning the number of terms retained in the series (2,4). 

The general notions concerning the increase in L as Re increases are 
of little help* as we shall see subsequently, To clarify this ques- 
tion for several values of Re, experimental calculations of the cur- 
rents were performed for various values of L from II to III; the re- 
sults of the calculations for Re * 1300 and 1500 in the form of the 
spectral component of the current velocity with respect to A for 
various values of L are 1 presented in Figures 5 and 6. 

The difficulty in choosing L, as is clear from Figures 5 and 6, 
is due to the fact that in critical situations (loss of stability 
of the basic current), the spectrum drops rapidly when L > 10, so 
that nothing is determined. Indeed, when Re » 1300, the spectra of 
the r-component of the velocity when L = 11, 21, and 31 are practi- 
cally identical, and the order of the next terms neglected does not 
exceed 1< of the maximum value. However, as the calculations with 
L - 39, 63 , 75, and 91 have shown, the actual spectrum has little in 
common with the spectrum for small L. The point is, as experiment 
shows, that near this value of Re the basic current becomes unstable 
the emerging secondary current has a substantially smaller charac- 
teristic linear dimension, and in the spectrum of the meridional 
component there appears a second maximum for large wave numbers, the 
prediction of whose occurrence on the basis of calculations with 
L < 30 turned out to be impossible. The spectrum of the r-component 
changed qualitatively with the increase in L, and the kinetic energy 
included in the radial component of the motions is now contained not 
only in the large scale motions of maximum size with & - 2, but also 
in the small scale motions with wave numbers & 3 * 30 to 46, i.e., for 
a sufficiently accurate description of the current when Re 1300* 
it Is necessary' to take & &*60« 
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Thus, to obtain a real picture of the current beyond the limit 
of stability (Re c < Re = 1300)-, it is necessary to Increase the 
number of terms of the series more than three-fold. At the same 
time, to describe the current when Re < 1300 it is sufficient, as 
calculations showed, to take L = 11, It must be emphasized that five 
terms of the series (2,4) give smooth, rapidly determined solutions 
of the equations even for very large Re (calculations performed up 
to Re = 3000) when, as is familiar from experiment, [5, 6], not a single 
stability limit is traversed and the basic current simply does not 
exist, and the ratio of the energy of several of the first neglected 
terms to the total energy does not exceed the order of several 
percent. 

As the spectrum of shows for Re = 1500 (Figure 6), the 
qualitative picture of the current is easily discerned when L - 6l; 
when L = 75* 91* HI* the spectral characteristics of the current 
practically coincide. 

As noted above, reorganization of the solution was observed /19 

when Re > 1300, and the current emerging here corresponded to the 
formation near the equator of four annular axisymmetric vortices 
(Figures 7a, 8a, 9a), This solution exists up to Re - 2200; when 
the number Re increases further, there occurs a- further reorganiza- 
tion of the current, as a result of which another secondary current 
emerges with eight annular vortices parallel to the equation (Figures 
9c, 9c), The solution with 8 vortices exists in our calculations 
in the interval 1750 < Re < 2600; when Re = 2700, a solution did not 
appear for a * steady-state regime, and an oscillating periodic regime 
was established with period approximately equal to the period of the 
single rotation of the sphere. The cause of such oscillations is 
under study . 

Henceforth for brevity in diagram 4, and from time to time in 
the text, we shall use the following designations for the current 
regimes; we shall designate the current with 2i steady-state 

annular vortices the regime i (.1 * 1, 2, *3, 4), and we shall call 

/ 
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the basic current the regime 0 ;» the transition from one regime to 
another will be indicated by arrows. 

When Re decreased from 1700, the number of annular vortices in 
regime 4 decreased to 6 (Figures 7b, 8b, 9b). The regime 3 thus 

t 

obtained exists U up ir to Re * 2700, whdre it — like regime 4 — 
becomes oscillatory. 

For decreasing Re when Re » 1550, regime 3 goes over into 
regime 2, and when Re decreases still more, it become regime 0 (for 
Re * 1300), by-passing regime 1. 

To obtain a solution with two annular vortices (regime 1) which, 
as was known from experiment, arise when Re > 1212 ± 25 during loss 
of stability by the basic current [6], a special disturbance is 
introduced into the basic current whose form is close to that of 
the desired current. In Figures 7c, 10, and 11, the function of 
the flow and the values of the angular velocity on the surface 
r = 1.049 are shown, as are also the spectral characteristics of /20 

the azimuthal and meridional motions for two different current 
regimes (regimes 0 and 1) for the same Re ** 1270, calculated for the 
same number of spherical functions L « 91- 

Regime 1 exists in the interval 1260 < Re < 1400; for larger 
values of Re, the current is reorganized into the regime with four 
annular vortices, previously described (regime 2). When Re « 1260, 
the regime is reorganized into the basic current. 

Besides -the flow functions and the angular velocity of rotation, 
spatial spectra of the total kinetic energy and the individual 
components of the velocity were obtained. In Figure 10 - 19, the 
graphs of the change in the flow functions and the angular velocity 
with respect to 6 are shown in the center of the annulus when r - 
1.049; also included are the Fourier coefficients, U«, W* of the 
the meridional and azimuthal component of the velocity for various 
regimes of the current when Re ** 1270, 1600, 1900, and 2000* In 
Figure 20, the flow lines and a photograph of regime 4 when Re * 1900 

i 

; 
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a.3?e shown. The spectra of the .total kinetic energy for all regimes 
for Re « 127Q* 1600, 1900, and 2200 are presented in Figures 21 - 22* 

An anlysis of the energy spectral show that if we were interested in 
the basis kinetic energy of the motions, then clearly it would be 
possible in many cases to limit ourselves to L ^ ^0. However, here 
the meridional motion, whose basic scale close to the equator is 
much less than the scale of the azimuthal motions, *ould be greatly 

distorted. 

i 

i 

; 3. THE HYDRODYNAMIC INTERPRETATION OF THE RESULTS 

OBTAINED . COMPARISON WITH EXPERIMENT 

The numerical solutions of the problem posed in the thick 
annulus with 6 - 1.33 and e * -1 have demonstrated the uniqueness 
of the solution of the axlsymmetric current in the interval of Re 
considered. 

The conformity with experimental and numerical results pre- /21 

viously obtained [1 - 3, 8 - 10], the meridional current is concen- 
trated in the lower latitudes (Figure 1) as Re increases, in the 
profile of the angular velocity there appears a point of inflection, 
and the spectrum of the kinetic energy shows that the basic part 
of the energy is contained in the large-scale motions. 

The graph of the torque M as a function of Re (Figure 2) has 
no breaks up to the value of Re ^ 410', but near Re ^ 54 there is a 
point of inflection. Up to Re - 130, the experimental and numerical 
values of M agree with very good accuracy within 0.5£* For Re - 
150, 200, and 250, the calculated values of M are given by circles in 
Figure 2; as is clear from the figure, the difference between the 
measurements and the calculated values is on the order of 2.5&* 

recorded when Re = 250, and It Is possible in this connection that 

i 

for this value of Re the accuracy of the current .calculations 

declines for L * 31. Analysis of the spectra. of the kinetic energy 

near the breaks noted by Munson £4] for Re 0 and Re c also discloses 

**1 2 

the ‘absence of any reorganisation In the current. 
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The results discussed together with the results of experiment 
[3$ 113 make it possible to state a conclusion concerning the 
stability of a Couette spherical current in thick annuli! in the 
interval for He which was considered* 

In the thin annulus with <5 * 0.11 and e * -1, the basic axisym- 

metric Couette current exists up to Re » 1300 according to the cal-* 

! 

dilations s which is somewhat larger thbn the first experimental value 
of the critical number Re = 1212 ± 25. For values of Re > 1260, 
only a special type of disturbance makes is possible to obtain numeri- 
cally a current regime with two annular vortices. This impossi- 
bility of obtaining regime 1 from the basic current by a simple in- 
crease in Re is connected, obviously, with the condition of symmetry 
of the current relative to the equator which we imposed. In an 
experiment [6], the appearance of regime 1 was connected with the /22 
momentary disruption of such symmetry at the limit of stability. 

The cause of the approximately 4$ divergence In the numerical ■ 
and experimental values of the first critical value of Re (Re c * 

1260 by calculation Re c « 1212 ± 25 by experiment) is still not 
clear to us. 

Just as in the case of experiment, calculations show that for 
supercritical values of Re there exist several different steady- 
state current regimes at one and the same point in space. On the 
diagram (Figure 4a), the regions of existence (relative to Re) 
corresponding to the regimes of the calculations are denoted by 
lines opposite the numerals 0, 1, 2*, 3, and 4. The arrows Indicate 
the transition from one regime to another as Re Increases without 
the introduction of additional disturbances, a dotted arrow denotes 
the passage from the basic current to regime 1 when disturbances 
of special type are Introduced into the current. The upper mark 
corresponds to the domain of existence of these regimes in the 
experiment [6]. The wavy form of the mark Indicates the appearance 
of sinusoidal disturbances in the annular vortices In the corre- 
sponding regimes. 
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Regimes 3 and as is clear from the diagram* exist fox* larger 
values of Re than in the experiment . Regime t fox 1 Re - 2150 and 
regime 3 for Re * S30Q are reorganised into essentially non- 
axisymmetric regimes with spirals [5, 6]V In the calculations*. 
both solutions for Re * 2700 become oscillations with a period 
approximately equal to the pex'iod of one revolution* For the 
present, it is difficult to say whether these oscillations are 
purely a calculational phenomenon* or the period axisymmetric re- 
gimes are solutions of the Navier-Stokes equations. If the latter 

i 

is true* then they must be unstable relative to tne non-axisymmetric 
disturbances, since.- in the experiment on spherical [6] and cylin- 
drical [173 annulii no periodic axisymmetric current was observed. 

Qualitatively, the forms of the currents for the regimes ob- /S3 
tained by calculation and by experiment compare favorably t also . 

Thus in regime iJ, the third annular vortex from the equator always 
occupies more space than the remaining vortices, in the experiment 
as well as in the calculations. Supercompression of the lines of 
flow and the formation of closed- vortices in the basic current 
(current of the "cat's eyes” type), obtained in the numerical calcula- 
tion (Figures 7* 8, 9, 2QB) were recorded in experiment [6] in the 
form of regions with large velocity gradients (the broad* dark belts 

pelow the annular vortices in Figure 20a). 

1 

*• 

CONCLUSIONS AND SUMMARY 

The results of the numerical investigation of nonlinear steady- 
state axisymmetric currents in a viscous incompressible liquid in a 
thick (6 « 1.33) and a thin (<S * 0,11) annulus when only the 
internal sphere rotates ( e = -1) in the indicated invervals of 
variation of Re make it possible to draw the following conclusions: 

1, In the interval 5 < Re < 250 in a thick annuius, thex^e was 

obtained (in agreement with experiment) a unique steady-state axi- 

symmetric solution — the basic Couette* spherical current. 

* 


* 
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2, The absence of breaks in the curves of the torques M(Ke) 

• (Figure 2) and the kinetic energy E( Re) (Figure 3)* and also the;, 
absence of qualitative changes in the spectrum for the kinetic energy 
of the current' as Re increases* make it possible to draw an infer- 
ence concerning the stability of the Gouette spherical current in a 
thick annulus relative to the steady-state axisymmetric disturbances 
in the Re interval considered* Taking into account the results of 
experiment til 3> this conclusion can be generalised to the stability 
of the Gouette spherical current in thick annuli relative to arbi- 
trary disturbances when Re <; 410, 

3, In a thin annulus, only those steady-state axisymmetric 
current regimes are obtained numerically which were observed experi- 

; s . ■ i 

mentall, i*e., the basic current and currents with two, four, six, /2t 
and eight annular vortices , The regions of existence of the numeri- 
cal solutions, as is clear from Figure) 4, can be distinguished from 
the regions obtained experimentally. The prolongation of certain 
current regimes into the region of larger Re values can be explained 
by the fact that the conditions of axisymmetry and symmetry relative , 
to the equator are imposed on the motions, and the current regimes 
observed in these regions are essentially non-axlsymmetric and non- 
symmetric with respect to the equator, 

4, The basic steady-state axisymmetric Gouette current in a 

thin annulus exists up to the value Re = 1300, which is somewhat 
larger than the experimental value of the first critical Reynolds 
number Re ~ 1212 + 25, above which the basic current turns out to 
be unstable, * 

5, In the supercritical region, the numerical as well as the 
physical experiment, gives a non-unique solution for the same Re 
values, 

6, In a thin annulus, the spectral characteristics of the 

velocity components disclose the following features of the secondary 
currents: a large part of the energy! of differential rotation 

during supercritical regimes is contained in maximum scale motions 
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(the wave numbers l 2 - q) , ahd the energy of the meridional 
motions have two approximately equi valued maxima for & £ 2 and 
■X 40> i.e-,,. large-scale and small-scale motions which are equally 
justified from the standpoint of energy, *‘Xt is interesting to note 
that as He increases within the very same regime, the predominant 
scale of the small-scale part of the periodic motions increases, 
I.e*, the wave number £ decreases. 


Thus, numerical modeling of the nonlinear Couette spherical 
current and the secondary currents which emerge for the supercritical 
values of He produced results which are in good agreement with 
experiment. As a supplement to a physical experiment, the calcula- 
tions make it possible to determine the total kinetic energy of 
currents, its spectral distribution, to study the fine structure of 
currents, and in disputed cases to resolve the question of the 
stability of currents. 


APPENDIX 


In the calculations, the time-related quantities for the 
(n + l) 3 ^ annulus were obtained after completion of the calculations 
for the annulus. At this time there were files for the quanti- 

ties U (fc), 1 4 WeX*-) In the storage of the computer by means 

of which, with the aid of formula (2.3), the values 


(*> 9 ii- 

means of these quantities, the nonlinear terms ( k ?fy> £) > 


were calculated. By 
h 


I 


1, 2, 3 were computed. 


To reduce the time to calculation 


cr 


at the beginning before undertaking the problem, three imbrices were 
calculated and the results stored In the memory of the computer. The 

elements of these matrices are (in the generally accepted 

notation) : 


& e °(cxy>9i 


v :'* 

> i ” ' > 

V 

* — v ■ 4 ^— rf 

.2. J 

j Q 4. 

. _ , L~i 

* * * 1* 5 

1 > J 2 . 


y (with respect to 


l-{ 
'?. x ’ 


;(Aoi) 


the inverse f the matrices (avl) were stored in the computer. 


After calculating the values of < 5 ^ % 9 the difference 

equations ( 2 , 8 ) and (2,9) were solved. Here, ? \/e y /“^'are found 

from equations 1, 2 and 4 with the help of a matrix trial run, and 
V\/l is obtained with the aid of the usual run. (We draw attention 
to the fact that no difficulties are encountered here because of the 
absence of boundary conditions for p.) 


The major part of machine time was spent in calculating J"? 
and <p>' . 

I 

Choice of the function f is quite arbitrary. We assume: 
i - V 3 4- Ct -ret 

"- V. 5 " (&+<>) ^ 
v, . />*■- A • A = — i-t- c. ) 

£ § r 0 / ^ > ■ v ' 

The computed points ¥>[*£•') on khe boundaries £ I and 
£ r: /-f-S* were thicker here than in the center of the annulus. In 
calculations relative to s the number of points along the 

radius v£T = 20 * and in calculations for » 0,11 , we .assumed 


The calculations for 5 » T, 33 Were carried out with the 
values of the calculated parameter t * 0.6; when 6 « 0.11, it was 
assumed that v * Q.2, j 

The time to determine the solution depended on the initial data 
and the value of Re, Here near the critical values of Re, at which 
the current of the given regime lost stability, the time to obtain 
the solution increased. Thus, e.g, , when Re « 1750 with initial 
values corresponding to Re « 1800, the fourth sign for torque M was 
determined after 400 steps with respect to t, and for Re - 1700 and 
initial values corresponding to Re - 1750, 1200 steps were required. 

Since we are interested in cases in which the solution of the 
steady-state equations Is not unique, we devoted special attention to 
the extent to which these equations are effectively satisfied. Let 
us consider the following figures: when Re * 1250- in the first of 

(O ^ /j 

the equations, ^ »Q,III75S7I *‘XG ; here the contribution of 

the energy terms - 0,9086*1956; that of the viscosity terms - 

-0.00*16*123005, and of the pressure gradient - -0,086222768. 

The processing and development of the results were carried out 
according to separate programs. Since the solutions are presented 
In the form of interpolational sums, we were able to calculate the 
values of the desired functions at any point 8, which made it pos- 
sible to obtain a detailed picture of the current In the equatorial 
region (Figure 7-9). 
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CAPTIONS TO THE FIGURES 


1. Flow functions of the meridional current and the lines of /3Q 
constant angular velocity in the annulus 6 = 1*33# Re * 200 


2, The curve of the torque M(Re) for the annulus 6 « 1*33 
Experiment and theory O 

m 9 At 3 ' z* &{/■ 

3* The curve of the kinetic energy E(Re) for the annulus 
5 “ 1.33 


4. Existence diagram for various regimes of currents (a) and 
the wave length of a vortex near the equator as a function 
of the number Re (B) In the thin annulus 6 = 0.11 

5* The spectrum U&(&) for various L when Re = 1300 in the 
annulus 6 « O.il 


6. The spectrum U^(£) for various L when Re * 1500 In the 
annulus <5 = 0.11 

7* The current functions ii>*10 2 and the lines of constant 
angular velocity in the annulus 6 » 0.11 

a) Re = 1600 regime^ 2 3 

b) Re = 1600 regime' 3* 

c) Re - 1270 regime 1 

8. The same as Figure 7 with Re - 1900 

a) regime 2 S 

b) regime 3# 

c) regime 4 

9. The same as Figure 7 with Re = 2200 

a) regime 2, 

b) regime 3# 

c) regime 4 


10. The current function and the angular velocity as functions 

of 6 when r - 1.049, and the spectra of and Wp(&) 

for the basic current when Re - 1270 

11. Ditto for regime 1 when Re = 1270 /31 

12. Ditto for regime 2 when Re = 1600 

13* Ditto for regime 3 when Re - 1600 

14. Ditto for regime 2 when Re = 1900 

15* Ditto for/ regime ‘ 3 when rJ - 1900 

16. Ditto f*or regime 4 when Re = 1900 
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Figure 17. 
Figure 18, 
Figure 19* 
Figure 20* 

Figure 21, 

Figure 22* 

Figure 22. 


Ditto for regime 2 :when Re 2200 

Ditto for regime 3 .when Re '= 2200 

Ditto for regime 4 when Re = 2200 

The stream lines if>-10 2 and a photograph of the current of 

region 4 when Re = 1900 and 6 ® 0,11 

t 

The spectra of the kinetic energy of the currents in the 
thin annulus 5 ■» 0.11 


Figure 23 
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